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Abstract. We study the Gross-Pitaevskii equation with a repulsive delta function poten- 
tial. We show that a high velocity incoming soliton is split into a transmitted component 
'/^ ' and a reflected component. The transmitted mass (L^ norm squared) is shown to be 

in good agreement with the quantum transmission rate of the delta function potential. 
p I ' We further show that the transmitted and reflected components resolve into solitons plus 

, dispersive radiation, and quantify the mass and phase of these solitons. 

, 1. Introduction 

We study the Gross-Pitaevskii equation (NLS) with a repulsive delta function potential 

■ J idtu + \dlu — q6o{x)u + = 



(N 
O 

o 



u{x, 0) = uo{x) 

As initial data we take a fast soliton approaching the impurity from the left: 

(1.2) Mo(x) = e™^sech(a; - xo) , > 1 , a;o < 0. 

^ . Because of the homogeneity of the problem this covers the case of the general soliton 
^ ', profile Asech.{Ax). The quantum transmission rate at velocity v is given by the square of 
the absolute value of the transmission coefficient, see ()2.2p below, 

(1-3) T,{v) = Uv)\' = . 

. For the soliton scattering the natural definition of the transmission rate is given by 
(1-4) T^,{v) = hm = ^ hm ||n(t)Uo \\h , 

provided that the limit exists. We expect that it does and that for fixed q/v, there is a 
a > such that 

(1.5) T^{v) = Tg{v) + Oiv""), as u ^ +00 . 

Based on the comparison with the linear case (see (j2.2H) below) and the numerical evidence 
we expect (jl.Sp with a = 2. 

Towards this heuristic claim we have 
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Figure 1. Numerical simulation of the case q = v = 3, xq = —10, at times 
t = 0.0, 2.7, 3.3, 4.0. Each frame is a plot of amplitude \u\ versus x. 



Theorem 1. Let 6 satisfy | < 5 < 1. Ifu{x,t) is the solution of ()1.H) with initial condition 
()1.2|) and Xq < —v^~^, then for fixed q/v, 



(1.6) 

uniformly for 



\u{x, dx 



x>0 



y2 _|_ q2 



0{v 



as V — > +CXO , 



We see that by taking b very close to 1, we obtain an asymptotic rate just shy of v 
More precisely, we show that there exists 
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Figure 2. A plot of the numerically obtained transmission T^{v) versus 
velocity v for five values oi a = q/v = 0.6, 0.8, 1.0, 1.2, 1.4. The dashed lines 
are the corresponding theoretical v — > +oo asymptotic values given by 1/(1 + 



diverging to +oo as 5 f 1 and q/v — > +oo, such that for fixed q/v, ii v > Vq, then 

1 
2 



\u{x, t)\'^dx 



< CV^~2\ 



x>o ^ + Q 

The constant c appearing here is independent of all parameters (g, v, and 6). 

We have conducted a numerical verification of Theorem^- see Fig. |21 It shows that the 
approximation given by p.6p is very good even for velocities as low as ~ 3, at least for 

0.6 <a = q/v < 1.4. 

A more elaborate numerical analysis will appear in our forthcoming paper jU]. 

Our second result shows that the scattered solution is given, on the same time scale, by 
a sum of a reflected and a transmitted soliton, and of a time decaying (radiating) term - 
see the fourth frame of Fig. ^ This is further supported by a forthcoming numerical study 
inj. In previous works in the physics literature (see for instance |2I) the resulting waves 
were only described as "soliton-like" . 
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Theorem 2. Under the hypothesis of Theorem^ and for 

— + 1 <t< (l-6)\ogv, 

V 

we have, as v ^ +oo, 

u{x,t)=UT{x,t)+UR{x,t) + OL^ (^{t-\xo\/v)-^^^^ +Oli{v'- 
(1-7) ut{x, t) = ^^'PT^M{A\-.^)ti2j^^ sech(AT(x - xo - tv)) , 

ur{x, t) = e^^fle-'""+'(^«-"')*/2^R sech{An{x + xo + tv)) , 
where At = {2\tg{v)\ - 1) + , Ar = {2\rg{v)\ - 1)+, and 

ifiT = a.Tgtg{v) + M\tgiv)\) + (1 - A|)|xo|/2t; , 



1 



(fR = aigTgiv) + ipo{\rq{v)\) + {1 - Aji)\xo\/2v, 



Here tg{v) and r^iv) are the transmission and reflection coefficients of the delta-potential 
(see fj2.2|) ). When 2\tq{v)\ = 1 or 2\rq{v)\ = 1 the first error term in (jl.Tj) is modified to 
OL^{{\og{t-\xo\/v))/{t-\xo\/v))-^). 



Here and later we use the standard notation 

a'' a > , 
a < 0, 



:i.9) 



This asymptotic description holds for v greater than some threshold depending on q/v 
and S, as in Theorem ^ The implicit constant in the 0^2 error term is entirely independent 
of all parameters {q, v, and 6), although the implicit constant in the Ol-^ error term depends 
upon q/v, or more precisely, the proximity of and to |. 

A comparison of the transmission and reflection coefficients p.3|) of the 5 potential, and 
of the soliton transmission and reflections coefficients ()1.8|) . is shown in Figure El 

Scattering of solitons by delta impurities is a natural model explored extensively in the 
physics literature - see for instance |21;[H1; and references given there. The heuristic insight 
that at high velocities "linear scattering" by the external potential should dominate the 
partition of mass is certainly present there. In the mathematical literature the dynamics 
of solitons in the presence of external potentials has been studied in high velocity or semi- 
classical limits following the work of Floer and Weinstein |B] , and Bronski and Jerrard PP - 
see |7] for recent results and a review of the subject. Roughly speaking, the soliton evolves 
according to the classical motion of a particle in the external potential. That is similar to 
the phenomena in other settings, such as the motion of the Landau- Ginzburg vortices. 

The possible novelty in and (jl.7|) lies in seeing quantum effects of the external 

potential strongly affecting soliton dynamics. As shown in Fig. |2l Theorem ^ gives a very 
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Figure 3. Comparison of linear and nonlinear scattering coefficients as 
functions of a ^= q/v. 

good approximation to the transmission rate already at low velocities. Fig. ^ shows time 
snapshots of the evolution of the soliton, and the last frame suggests the soliton resolution 
fll.7|) . We should stress that the asymptotic solitons are resolved at a much larger time - 
see p. 

The proof of the two theorems, given below in ^BHU proceeds by approximating the 
solution during the "interaction phase" (the interval of time during which the solution sig- 
nificantly interacts with the delta potential at the origin) by the corresponding linear flow. 
This approximation is achieved, uniformly in g, by means of Strichartz estimates estab- 
lished in ^ The use of the Strichartz estimates as an approximation device, as opposed 
to say energy estimates, is critical since the estimates obtained depend only upon the 
norm of the solution, which is conserved and independent of v. Thus, v functions as an 
asymptotic parameter; larger v means a shorter interaction phase and a better approxima- 
tion of the solution by the linear flow. Theorem 2 combines this analysis with the inverse 
scattering method. The delta potential splits the incoming soliton into two waves which 
become single solitons. 
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2. Scattering by a delta function 

Here we present some basic facts about scattering by a delta-function potential on the 
real hne. Let g > and put 

2rfa;2 

We define special solutions, e±{x, A), to {H^ — /2)e± = 0, using notation given in (jl.yj) : 
(2.1) e±(x, /) = tg{l)e^'^^xl + (e^^'^ + rg(/)e^*'^)x° 



^1 = -n— + 9^o(^)- 



where tq and Vq are the the transmission and reflection coefficients: 

(2.2) «.W = a^. '-'(^'^a^- 

They satisfy two equations, one standard (unitarity) and one due to the special structure 
of the potential: 

(2.3) \tq{X)\' + \rq{l)\' = 1 , tq{X) = 1 + rq{X) . 

We use the representation of the propagator in terms of the generalized eigenf unctions- 
see for instance the notes [15] covering scattering by compactly supported potentials. The 
resolvent 

RqiX) (Hq - X' /2)-\ 

has kernel given by 

Rqil){.x,y) = — ^ {e+{xj)e^{y,l){x - y)l + e+{yj)e^{x,l){x - yf_) . 

lltq{l) 

This gives an explicit formula for the spectral projection, and hence the Schwartz kernel of 
the propagator: 

(2.4) eM-^tHq) = — e-*^'/2 (^e+(x, /)e+(y, /) + e_(x, /)e„(y, /)) dl . 

The propagator for Hq is described in the following 

Lemma 2.1. Suppose that ip E and that suppv? C (— oo,0]. Then 

(2 5) ^^'Pi-'^tHq)^{x) = 

exp(-itifo)(v5 * Tq){x)x'l + {exp{-itHo)(p{x) + exp(-iti7o)(v^ * Pq){-x))x^_ , 

where 

(2.6) pq{x) = -qexp{qx)x°_ , Tg{x) = 6q{x) + pq{x) . 
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Proof. All we need to do is to combine ()2.H) and ()2.4p . Using the support property of (f we 
compute, 



(p{y)e+{y, \)dy = rq{-X)ip{-X) + 0{X) 



so that 

f oo 



(exp(-ztff,)^)U>o = ^ / e-''^'/' (t^(A)e^^-(r,(-A)^(-A) +^(A)) + 







(r,(A)e*^^' + e-*"^)t,(-A)<^(-A)) dX 
= — [ e-''^'/Hg{X)0{X)e'^'^dX 

= exp{-itHo){Tq * (p){x) , Tq{X) = tg{X) , 

where we used the fact that rg(— A)tq(A) + rq{X)tq{—X) = 0. 
Similarly, using rg(— A)rg(A) +tg(— A)tg(A) = 1, we have 

1 f°° 

iexp{-ttHg)if)\,^o = — / e-^*^'/2 ^0(^x)e'Xx ^ r,{X)0iX)e-'^^) dX 

27r Jo 

= exp{-itHo)(p{x) + exp{-itHo){(p * Pq){-x) , pq{X) = rg(A) . 
A simple computation gives ()2.(ij) concluding the proof. □ 

We have two simple applications of Lemma I^ITI the Strichartz estimates fProposition l2.2|) 
and the asymptotics of the linear flow exp{—itHg) as v ^ +oo fProposition I2.3|) . We start 
with the Strichartz estimate, which will be used several times in the various approximation 
arguments of ^ Since it is particularly simple in our setting, we give a complete proof 
(see for references and the general version). 

Proposition 2.2. Suppose q > and 

(2.7) idtu{x,t) + ^dlu{x,t) - q5o{x)u{x,t) = f{x,t) , u{x,0) = (p{x) . 

Let the indices p,r, p, r satisfy 

, , 2 11 2 15 

2.8) 2<p,r<cx), l<p,r<2, - + - = -, ^ + ^ = 77 

p r I p r 2 

and fix a time T > 0. Then 

(2.9) lkllLf„,,Lj <c||^|U2 + c||/||^, , 

[0,1 \ ^ [0,T] a: 

The constant c is independent of q and T. Moreover, in (I2.7|l . we can take f{x,t) = 

-p jf with \\q\\ - 



g{t)So{x) and, on the right-hand side of ()2.9j) . replace \\f\\i^p ^j^f with \\g\^_^^ 



[0,T] 
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Proof. We put Uq{t) =^ exp{—itHq), so that Uq{t) is a unitary group on L^(M). For cp e 
L^(M) we have, using Lemma f2.H 

\\Ug{t)ip\\L^ <J2\\Ugit){vxl)\\L^ 
± 

< E ro(t)IUi^L^(ll(^4) * Tqh^ + 11(^4) * p.lUO 



(2.10) 



< 



< 



3 



:(l+2||p,|Ul)||(^|Ul 



By the Riesz-Thorin interpolation theorem (see for instance fW\ Theorem 7.1.12]) we have 



(2.11) \\U(t)\\^.,^^.<C\t\-^^'-^), l<r'<2, ^ + ^ = 1 

The estimate ()2.9|) with / = reads 



\\um\L^,Lr<C\\g\\LHM) 



which by duahty is equivalent to 



(2.12) 



U{-s)F{s)ds 



< C\\F\ 



L2( 



tp' Tr' • 



The two equivalent estimates together give f ()2.12p is applied with p', r' replaced by p, r 
it is easily checked that ()2.8p still holds) 



[ U{t-s)F{s)ds 
Jr 



< \\F\ 



Putting F{s) = ]l[o,t](s)/(s, x) we obtain ()2.9|) for uq = 0. Hence it suffices to prove ()2.12|) . 
Put 

TF{x)= [ U{-s)F{s,x)ds, 
Jr 

and note that T*g{s,x) := U{s)g{x). The estimate ()2.12|) is equivalent to 
which is the same as 



(2.13) 



( [ (Git), U{t-s)F{s)) dtds 



<C\\G\\,,,AF\\,,,,. 
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To obtain ()2.13p from (j2.11|) we apply the Hardy-Littlewood-Sobolev inequality which 

says that if Ka{t) = and 1 < a < oo then 

111 

||i^a * i^|U"(R) < , - = - , l<(3<a, 

cx jO a 

see for instance fOl Theorem 4.5.3]. We apply it with 

a 2 V 

which is the admissibility condition (j2.8|) . □ 



a = p , (3 = p' 



We now turn to the large velocity asymptotics of the linear flow exp(—itHq). 

Proposition 2.3. Let 9 G C°°(]R) he hounded , together will all of its derivatives. Let 
if G V > 0, and suppose supp[^(«)v9(« — Xq)] C (— cx3,0]. Then for 2\xq\ / v <t<l, 

(2.14) e-'*^« [e'"V(x - xq)] = t(t;)e-'*^° [e"V(x - xq)] 

+ r(t;)e-^*^°[e-"X-x - xo)] 
+ e{x,t) 



where, for any k > 0, 



|l2 < -\\d^[6ix)(f{x - Xo)]\\l2 



V 

(tvY 

+ mi-e{x))if{x-xo)\\Li 



+ 7Z:TllK^)V(a;)||Hfc 



In ^ Proposition 12.31 will be applied with 6{x) a smooth cutoff to x < 0, and v^(x) 
sechx with xq = —v^~^ <^ 0. 

Before proving Proposition 12.31 we need the following 

Lemma 2.4. Let tjj G 5(M) with supp^ C (— oo,0]. Then 

(2.15) e-**-^'[e"'Xx)](x) = e-^*^°[e"Xx)](x) x° 

+ t(i;)e-**^°[e*"Xx)](x) x° 
+ r(t;)e-**^°[e-"X-x)](x) x° 
+ e(x, t) 



|e(x,t)||i2 < -||(9^^||l2 



uniformly in t. 
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Proof of Lemma \2.4\ By ()2.5|) with (f{x) = e^^'"ip{x), 

e{x,t)= [e-*^o(y,*(r-t(^)5o))(x)]x° 

+ [e-'*^«(v^ * (p - r{v)6o))i-x)] x° 

and thus it suffices to show 

(2.16) ||e-^*^°(^ * (r - t{v)5o)){x)U2 < -Wd^^Ll 



and 



|e-^*^°(^*(p-r(t;)5o))(x)|U. <-||a. 



V 



The proofs of these two estimates are similar, so we only carry out the proof of ()2.16|) . 
By unitarity of e~^^^° and the Plancherel's identity, 



(2.17) ||e^*^°[<^ * (r - t(t;)5o)](x) |U| = U{\ - v){t{\) - t{v) 

Since 

—iq{X — v) 



t(A) - tiv) 



{iX — q){iv — q) 

we have \t{X) — t{v)\ < {X — v)/v. Using this to estimate the right-hand side of (j2.17|) and 
applying Plancherel's identity again yields ()2.16p . □ 

Proof of Proposition 1^31 Apply ()2.15|) to ip^x) = 6{x)(p{x — Xq) to obtain 

(2. 18) e"**^' [e*"V(x - xo)] (x) = e"''*^° [e*"V(a; - xq)] (x) x° 

+ t{v)e-''"° [e*"V(x - xo)] (x) x° 
+ r{v)e-''^° [e-"'V(-x - xq)] (x) x° 
+ ei(x, t) + e2(x,t) 

where ei(x, t) is as in Lemma ITU and (putting /(x) = e*'^*'(l — 6'(x))(/9(x — Xq)) 

62 (x, t)= + e-''^^f{x) - e-**^V(a;) 

- t(w)e-**^"/(a;) - r(w)e-^*^° [/(-x)](x) x° 

By Lemma f2. 41 

||ei(x,t)||L2 < ^||a^[6l(x)v3(x - xo)]||l| 
uniformly for all t, and by unitarity of the linear flows, 

l|e2(a;,t)|U2 < 4||(1 - ^(x))^(x - xo)|U2 
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also uniformly in all t. Now restrict to the time interval 2\xq\/v < t < 1. By ()2.18|) . it 
remains to show that 



(2.19) 



(tv) 



Ck 



e-*^«[e-X-x-Xo)](x)|L.^^ < 



The second of these is in fact equivalent to the first, since for any function g{x), 

e~^*'^'>[g{-x)]{x) = e-^'''^[g{x)]{-x) . 
Now we establish ()2.19|) . Since 

[e^-V(. - XoWiX) = e-"°(^-^)^(A - v) , 

1 



e-**^«[e^^Xx-Xo)](x) = — / e— e — -e 



271 

e 



2tt I 



By k applications of integration by parts in A, 



^X — Xq — tv ^ 

Since 2\xo\/v < t, we have —xq — tv < and thus \x — xq — tv\ > \ — xq — tv\ > tv/2 for 
X < 0. Hence 

|e-'*^o[e'"V(x-xo)](x)|2rfx 



(2.20) 



< 



Ck 



(tv)'' 

Ck 

(tv)'' 



from which the result follows by applying the Leibniz product rule and the Plancherel 
identity once again (and using that t <1). □ 



Remark. Suppose that 

m(x, t) = e-''^" [e^^Xa;)] , V e '5(M) , supp C (-oo, 0) , 
Then for t ^ 1 and as f ^ +oo, 

,2 / I 



L2 



1. 



(2.21) 



|M(x,t)| dx 



y2 _|_ q2 



V 
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In fact using (j2.5|) and an estimate similar to (j2.20|) we see that for t > 1 

/ |M(x,t)|2rfx= ||e-^*^o((e^*»*r,)x^f = ||(e*«»*r,||2 + 0(t;— ) 
Jo 

Ztt 

An expansion in powers oi {X — v)/v gives ()2.2H) . 

3. SOLITON SCATTERING 

In this section, we prove Theorem^ We recall the notation for operators from Sect|21 
and introduce short hand notation for the nonlinear flows: 

• Hq = —\d1. The flow e~^^^° is termed the "free linear flow" 

• Hq = —\d1 + q5Q{x). The flow e"**-'^'' is termed the "perturbed linear flow" 

• NLSg(t)<y9, termed the "perturbed nonlinear flow" is the evolution of initial data 

according to the equation idtu + \d1u — q6o{x)u + = 

• NLSo(t)</3, termed the "free nonlinear flow" is the evolution of initial data (p{x) 
according to the equation idth + ^d^h + \h\'^h = 

From Sect^we recall the form of the initial condition: Uq{x) = e*^''sech(a; — xq), v ^ 1, 
xq < —v^~^, I < 5 < 1, and we put u{x,t) = NLSq(t)uo(x). 

We begin by outlining the scheme, and will then supply the details. The O notation 
always means difference, uniformly on the time interval specified, and up to a multi- 
plicative factor that is independent of q, v, and 6 (any such dependence will be exhibited 
explicitly). 

Phase 1 (Pre-interaction). Consider < t < ti, where ti = \xo\/v — v^^ so that 
Xq + vti = —v^"^. The soliton has not yet encountered the delta obstacle and propagates 
according to the free nonlinear flow 

(3.1) u{x, t) = e-^*^'/2git/2gix«ggpj^^^ ^ 0{qe-^'''), < t < ti 

The analysis here is valid provided v is greater than some absolute threshold (independent 
of q, V, or 6). But if we further require that v be sufficiently large so that v~^^'^e'"^ ^ > a 
(recall a = q/v), then qe'""^ < < v~^l'^. This is the error that arises in the main 

argument of Phase 2 below. 

Phase 2 (Interaction). Let ^2 = + 2t'~'' and consider t\<t < ti- The incident soliton, 
beginning at position — f encounters the delta obstacle and splits into a transmitted 
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component and a reflected component, which by time t = t2, are concentrated at positions 
v^~^ and —v^~^, respectively. More precisely, at the conclusion of this phase (at t = ^2), 

(3.2) u{x, = t{v)e-''^'''/^e''^/^e'''''sech{x - Xq - vt2) 

+ r{v)e-''^'''l^e''^'^e-'''\ech{x + xq + ^^2) 

This is the most interesting phase of the argument, which proceeds by using the following 
three observations 

• The perturbed nonlinear flow is approximated by the perturbed linear flow for 

ti<t< t2. 

• The perturbed linear flow is split as the sum of a transmitted component and a 
reflected component, each expressed in terms of the free linear flow of soliton-like 
waveforms. 

• The free linear flow is approximated by the free nonlinear flow on ti < t < t2- Thus, 
the soliton-like form of the transmitted and reflected components obtained above is 
preserved. 

The brevity of the time interval [^1,^2] is critical to the argument, and validates the ap- 
proximation of linear flows by nonlinear flows. 

Phase 3 (Post-interaction). Let ts = ^2 + (1 — <5) logv, and consider [^2,^3]- The trans- 
mitted and reflected waves essentially do not encounter the delta potential and propagate 
according to the free nonlinear flow, 

(3.3) u{x, t) = e~''*"'/2e**^/2e""NLSo(t - ^2) [t{v)sech{x)] {x - xq - tv) 

+ e-^*'''/2g*i2/2g-ix,;^LSo(t - t2)[r{v)sech{x)]{x + Xq + tv) 

This is proved by a perturbative argument that enables us to evolve forward a time (1 — 
5) \ogv at the expense of enlarging the error by a multiplicative factor of e^^"''^'"^^ = v^~^ . 
The error thus goes from tr^^'^ at t = ^2 to t;-^~2<^ at t = ^3. 
Now we turn to the details. 

3.1. Phase 1. Let Mi(x, t) = NLSo(t)Mo(a;) and u{x,t) = NLSg(t)uo{x). Let w = u — Ui. 
Recall that ti = \xq\/v — so that Xq + vti = —v^~^. Note that 

We will need the following perturbation lemma. 

Lemma 3.1. Ifta < h, U - ta < ci, and \\w{-,ta)\\L^ + q\\ui{0,t)\\L?° , < 1, then 

IklU- ^ M < C2(|k(-,ta)||L2 +g|kl(0,t)|U- ^ ) 
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where the constants C\ and C2 depend only on constants appearing in the Strichartz estimates 
and are, in particular, independent of q and v. 

Proof, w solves 

idfW + d^w — q6o{x)w = — \w + Ui\'^{w + Ui) + + q6o{x)ui 

= — \w\'^w — (2ni|u;|^ + Uiw'^) — (2|ni|^u; + ulw) +q5o{x)ui 



-V ' ^ ^ 



cubic quadratic linear 

From this equation, w is estimated using Proposition 12.21 For the cubic nonhnear term we 
take p = r = 6/5 and estimate by Holder as 

IlkPw^lU/s r6/s < (tb-taY^'^WwWle || W || , ,L2 

Since complex conjugates becomes irrelevant in the estimates, both quadratic terms are 
treated identically. In Proposition 12.21 we take p = r = 6/5 and estimate by Holder as 

\\uiW^\\ e/5 ^e/s <{th-taY^^\\w\\le rslkllUff , 



< V2{t!, - taY 



/2| 



W\ 



For the linear terms (both of the form ulw), we take p = r = 6/5 in Proposition 12.21 and 
estimate as 

\\ulw\\ 6/5 6/5 < (tfe - taY^^\\w\\Lf reJUlllLl ^ ^L^JMl^ , ,Ll 



2/3 I 

The delta term is estimated by the concluding sentence of Proposition 12.21 as 

g||w(0,t)||^4/3 <q{tb-taY/^uiO,t)\\L^ 

Since tb — ta < 1, collecting the above estimates we have (taking \\w\\x = \\w\\l°° + 

\\w\\x < c\\w{-,ta)\\Ll + C{tb - taY^^i\\w\\x + \\w\\l + \\wfx) + CqMO,t)\\L^ 

Provided [t^ — taY^'^ ^ V (2c) above, the linear term on the right can be absorbed by the 
left as 

\\w\\x < 2c\\w{-,ta)\\Li + - 4)'/'(lkllx + Ik'llx) + 2cg||n(0,t)|Ufx> 

Continuity of ||tt^||x(4) as a function of t^, shows that provided 2c{tt — tay^'^{4:c\\w{-,ta)\\L^ + 
4cg|lwi(0, t)||L^ ) < 1/2, the above estimate impHes 

\\w\\x<^c\\w{;ta)\\L^^+4cq\\u{0,t)\\L^^^^^^ 



<2(t,-4)2/3||w;||^a g 
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concluding the proof. □ 

Now we proceed to apply Lemma 13.11 The constants ci and C2 will, for convenience 
of exposition, be taken to be ci = 1 and C2 = 2. Let k > he the integer such that 
k <ti < k + 1. (Note that /c = if the soliton starts within a distance v of the origin, i.e. 
—V — v^~^ ^ Xq < —v^~^, and the inductive analysis below is skipped.) Apply Lemma f3. II 
with ta = 0, tb = I to obtain (since w{-, 0) = 0) 

IklU-^jLi < 2g||ui(0,t)||L^^j < 2gsech(xo + tO 

Apply Lemma f3. II again with = 1, = 2 to obtain 

IklU^^Li < 2(||w;(-, + q\MO,t)\\L^J 

< 2^gsech(xo + v) + 2"'^gsech(xo + 2v) 

We continue inductively up to step k, and then collect all k estimates to obtain the following 
bound on the time interval [0, k] 

k 

||'"^||l-,jL2 < 2g^2''-^sech(xo + jt;) 
i=i 

The estimate secha < 2e~'"' reduces matters to bounding 

k-l 
j=0 

and, after summing the geometric series, we obtain 

\\w\\l^ l2< c2'=e"»+^'^ < cge"o+*=" 

" "-^[o.fe]-^^: — 2~^e'" — 1 ~ 

where the last inequality requires 2~^e'" > 2. Finally, applying Lemma f3. II on [/c, ti]. 



xo+kv 

As a consequence, ()3.H) follows. 



^Wl^^ ji? ^ ciqe'^^ + gsech(a;o + hv)) < cqe~ 



3.2. Phase 2. We shall need a lemma stating that the free nonlinear flow is approximated 
by the free linear flow, and that the perturbed nonlinear flow is approximated by the 
perturbed linear flow. Both estimates are consequences of the corresponding Strichartz 
estimates (Proposition I2.2|l . Crucially, the hypotheses and estimates of this lemma depend 
only on the norm of the initial data (p. Below, fl3.5j) is applied with (f{x) = u{x, ti), and 
||'w(3;, ii)||L2 = ||mo||l2 is independent of v; thus v does not enter adversely into the analysis. 

Lemma 3.2. If Lp ^ L"^ and < t{, such that t^ < Ci\\ip\\j^2, th 



en 



(3.4) ||NLSo(t)<^ - e-^'*^VlUfs=,,^jLi < c^tl^'Mh 



16 



J. HOLMER, J. MARZUOLA, AND M. ZWORSKI 



(3.5) ||NLS,(t)(^ - e-^*^vllLfj^,^,Li < C2CMh 

where c\ and C2 depend only on constants appearing in the Strichartz estimates. In partic- 
ular, they are independent of q. 

Proof. Estimate (|3.4|) is in fact a special case of (|3.5p obtained by taking g = 0. Let 
h{t) = NLSg(t)(^ so that 

idth + Idlh - q6o{x)h + \h\^h = 
with h{x,0) = f{x). Let us define 

with the natural norm, || • \\x- We apply Proposition 12.21 with, in the notation of that 
proposition, u{t) = h{t) — e~'^'^^''ip, f = —\h\'^h, p = r = 6, p = r = 6/5, and then again 
with p = oo, r = 2, p = r = 6/5, to obtain 

\\h{t) - e-^'^'^ipllx < C\\\h\''h\\^e/, ^6/5. 

The generalized Holder inequality, 

\\hih2h3\\p < II /ii II II II 92 II /is II 93 . ^ = :r + :r + ;r ' 

P Ql 12 (l3 

applied with hj = h, p = 6/5 and qi = q2 = 6, qz = 2, gives 

||Mt)-e-^*^'^|U<C||/.||ia^^^^^^e||/z|L.^^^^^^ 
(3-6) <Ctl/'\\h\\l, rML-^,Ll 



<Ctl/"\\h\\'x 



Another application of the homogeneous Strichartz estimate shows that 

||e-**^Vllx<C||(/.|U2, 

and consequently. 



x<c|iv^iU2+ct;/'ii/ii|3^ 

By continuity of ||/i|U(t6) in tb, if ctl^'^{2c\\(p\\L2y < 1/2, 

\\h\\x < 2c||v2||l2 

Substituting into (|3.6p yields the result. □ 
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Now we proceed to apply Lemma EUl Set ^2 = ^1 + St;"*^, and apply (|3.5p on [^1,^2] to 
obtain 

= NLS,(t-ti)[M(-,ti)] 

= e-'(*-*^)^'[u(-,ti)]+C(i;-''/2) 

By combining this with (jH.lll . 

(3.7) u{-,t) = e-'*^^''/2e'*i/2e-*(*-*i)^'[e'""sech(a: - xq - hv)] + 0{v^^/^) 

By Proposition 12.31 with 6'(x) = 1 for a; < — 1 and 9{x) = for x > 0, ip{x) = sech(x), 
and Xq replaced by xq + tiv, 

e-'(*^-*i)^'[e'"''sech(x - xq - vti)]{x) 

= t(t;)e-*(*^-*^)^o[e'"^ech(x - xq - vti)]{x) 
+ r(i;)e-^(*^-*^)^"[e-"'^sech(x + Xq + vti)]{x) 
+ 0{v-') 

By combining (jSHj), and (Q, 

u{.,t) = t(?;)e-'*^"'/2e'*i/2NLSo(t2 - ti)[e'"''sech(x - xq - 

+ r(i;)e-**i^'/2giti/2^LSo(t2 - ti)[e-^^^sech(x + + 

By noting that 

NLSo(t2 - ti)[e^^''sech(x - xq - hv)] 

and 

NLSo(t2 - ti)[e-*^^sech(x + Xq + hv)] 

= g-^(t2-tl)^'V2g^(^2-^l)/2g~^x.gg^J^^^ ^ ^ ^^^^^ 

we obtain ()3.2p . 

3.3. Phase 3. Let ^3 = ^2 + (1 — logv. Label 

utr{x,t) = e-**'''/2g*t2/2g*x^NLSo(t - t2)[t{v)sech{x)]{x -xo- tv) 
for the transmitted (right-traveling) component and 

Uref{x,t) = e"'*"'/2e**^/2g-ix.^LSo(t - t2)[r(^;)sech(x)](x + xq + tv) 
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for the reflected (left-traveling) component. By Appendix ^ for each /c G N, there is a 
constant c{k) > and an exponent cT{k) > such that 

c(k)(\os:vY^''^ 

(3.9) WuUxMlI^, + IMxMli^, + K{o,t)\ + Mo,t)\ < ^ '^^JJ^ 

uniformly on the time interval [^25^3]- We shall need the following perturbation lemma, 
again a consequence of the Strichartz estimates. 

Lemma 3.3. Let w = u — Utr — Wref- If ta < h, tb — ta < Ci, and 

c(fc)(g)(logtO-W ^ 

+ < 1, 

then 

II II < (w ( , \\\ , c(fc)(g)(logt;)-W 
ll^ll^l?...,^^ ^ ^2 (^lk(-,ta)|Ui + ^^^(Tri^ 

The constants Ci, C2 depend only on constants appearing in the Strichartz estimates and 
are in particular independent of q and v. 

Proof. We write the equation satisfied by w: 
idtw + \dl.w — q6o{x)w 

= - \W + Utr + MrcfP(w + Mtr + Mrcf) + IWtrpMtr + IWrcfpMref 

+ q6o{x)utr - q5o{x)u,ci 

= - \w\'^W - (2(Mtr + Mref)k|^ + (Wtr + Mref)^^) " (2|Mtr + Mrcfl^W + {Utr + U^cffw) 
- (Wtr^rcf + 2Urcl\Utr\'^ + u'^^fitr + 2Mtr l^rcf + 9^0 (^)Wtr - ^^0 (^)Wref 
tr— rcf interaction tr— delta ref— delta 

w is estimated using Proposition 12 . 21 and Proposition 12.21 The cubic, quadratic, and linear 
in w terms on the first line are estimated exactly as was done in the proof of Lemma 13.11 
For the "tr — ref interaction terms" (taking Mrcf|wtrP as a representative example), we apply 
Proposition 12.21 with p = A/3, r = 1 and estimate as 

(3.10) IIWreflMtrHl r4/3 ,i < c(t;, - 4)^^"^ || Mtr || Lff , ,L2 ||MrefMtr||Lff , ,L2 

||'"tr||L2 = "\/2|i('i^)| by mass conservation for the free nonlinear flow, and 

llMtrMreflU- ^^^jL2 < ||MtrMrcf _^^,l2^„ + || Mtr^ref || L^^^ 
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Now 

lkref|U-Li° = ||NLSo(t)[r(w)sech](x)||Loc^^ 

< ||NLSo(t)[r(t;)sech](x)||^/^^2||5,NLSo(t)[r(t;)sech](x)||^/^^2 

< c 

by mass and energy conservation of the free nonlinear flow. Similarly, ||iitr||Lf>Lg° ^ c. By 
this and ()3.9|1 . the above yields 

c{k) [log vY^''^ 

WUtrUrefh^^^^^^Ll < 

Thus, by (jSHHD, 

r-^in II I |2|| ^cimogvr^ 

and similarly for all other "tr — ref interaction" terms. Now we address the "tr— delta" and 
"ref— delta" terms (working with q6o{ the representative of both). By Proposition 

12. 2| we estimate as 

g||Wtr(0,t)||^4/3 <cih-tam\uUO,t)\\L^ 

By (El), 

(3.12) g||n,,(0,t) 11^4/3 <c{k)q{U-Uf' ^ I \ 

Collecting ()3.1ip . ()3.12|) . and the estimates for cubic, quadratic, and linear terms in w (as 
exposed in Lemma ITTl) . we have, with \\w\\x = l2 + llwlLe ^6 

\H\x < C\\W{; QU. + C(4 - Q'/'{\\w\\x + Ml + 11.^11^) + '^^^^^11!!^^'^^'^ 

If c(tb — ta)^^^ < 1 5 then the first-order w-term on the right side can be absorbed by the 
left, giving 

By continuity of in tb, if 

we have 

n n . n / xn 4c( fc) (g) (log w)<^(*^) 

Ikllx < 4c\\w{;Q\\l^ + ^ ' 

completing the proof. □ 
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Assume that a = q/v has been fixed. Choose k = k{5) large so that k{l — 5) > 3. Then 
the coefficient appearing in Lemma 13.31 is bounded by 



c(fc)(g)(logt;)"W ^ /q\ (logi;)'^W 



Now take v sufficiently large in terms of {q/v) and k (thus in terms of 5) so that the above 
is bounded by v~^. 

Now we implement Lemma 13.31 For convenience of exposition, we take Ci = 1, C2 = 2. 
Let £ be the integer such that £ < (1 — 5)logf < £ + 1. We then apply Lemma 
successively on the intervals [^2, ^2 + 1], • • • , [^2 + ^ — 1, ^2 + -^l as follows. Applying Lemma 
l3.3l on [^2,^2 + 1]; "we obtain 

\\w{-ML-^,^,,^Ll<m<.t2)\\Ll+V-') 

Applying Lemma f3. 31 on [^2 + 1,^2 + 2] and combining with the above estimate, 

^ n<M)\\Li + (2^ + 2)tri 
Continuing up to the £-th step and then collecting all of the above estimates, 

\\w{-Ml^^,^,lI < 2%w{-M)\\lI + (2' + ■ ■ ■ + 2)^-^ 
Since ||w(-,t2)IU| < v^^''^ and 2^ < v'^~\ 

(3.13) \\w{-Ml^^^,^,li<cv^-'^' 
thus proving ()3.3|) . 

Now we complete the proof of the main theorem and obtain ()L6|) . By ()3.13|) and ()3.9|) . 

(3.14) ||M(-,t) - u,,{-MlI^, < lk(->^)llL^>o + ll«-f(-'^)IU^>o ^ c^'"'' 

Since ||Mtr(-, i)||L2 = t{v), ()3.9|) implies ||Mtr(-, i)!^^^^ = ^("^0 + C(f "'^~^'^), which combined 
with ()3.14j) gives p.fip and proves Theorem ^ 

4. Resolution of outgoing waves 

In this section, we prove Theorem |21 We note that the proof of Theorem presented in 
§ni in fact provided a more complete long-time description of the solution: 

(4.1) u{x,t) = e-'*"'/2e'*^/2e'"''NLSo(t - t2)[t(t^)sech](x - xq - tv) 

+ e-'*"'/2e'*^/2e-'"''NLSo(t - t2)[r(i;)sech](x + xo + tv) 

+ Ol2{v^-I^), ^-^ + v~^ <t<c{l-6)logv 

where t{v), r{v) are defined in ()2.2p and NLSo(t)v2 denotes the solution to the NLS equation 
idth+ ^dlh+ \h\'^h = (without potential) and initial data h{x, 0) = (f{x). It thus suffices 
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to obtain the resolution of NLSo(t — t2)[t{v)sech.] and NLSo(t — t2)[r{v)sech] into solitons 
plus radiation decaying in L^. By the phase invariance of the free nonlinear flow 

NLSo(t - t2)[t{v)sech] = ^^NLSo(t - t2)[|t(^^)|sech] 

and similarly for NLSo(t — ^2) ['"(^')sech]. Since < k(^OI — ^-Pply asymptotics 

flB.l|) proved of Appendix IE) using the inverse scattering method. When \t{v)\ or \r{v)\ 
is equal to 1/2 we use the result of ^21 recalled in ()B.2|) . The result obtained by these 
substitutions differs from that stated in Theorem |21 by a factor of 

(4.2) expfz — 2^'^~^ 

for UT{x,t), owing to the fact that t2 = \xo\/v + v~^. But ()4.2j) differs from 1 by ~ v~^, 
and thus omitting it only introduces a discrepancy of in both and L^. There is a 
similar inconsequential disparity in the U{i{x,t) part. 

Appendix A. Spatial localization of the free nonlinear propagation 
Let ip E S and 

(A.l) 



idth + \dih+\hYh = Q, 
h{x, 0) = (p{x). 



Notational conventions. We denote by d hereafter. The x and t dependence of h{x, t) 
will be routinely dropped. The constants c{k), a{k), and the polynomials g-yit) that appear 
below may change (enlarge) from one line to the next without comment. The constants 
c{k) depend on the fixed function ip E S. 

The solution h satisfies conservation of mass and conservation of energy, which means 
that the integrals 



Eo= / \h\^dx, E2 = - ! {\dh\^ - \h\^)dx , 
Jr Jr 



are independent of time t. Since ||/i|||oo < ||/i||L2|l'9^|U2, we have ||/i|||4 < ||/i|||2 ||c?/i||L2 
and it follows from the E2 and Eq conservation that ||9/i||l2 < c, where c depends on Hv^Hlz 
and ||9v9||l2. 

In fact, there are an infinite number of conserved integrals, Ek, with integrands defined 
inductively as follows 



(A.2) 



^ ^ ii+i2=fc-i 



see 113 §8] for a proof of this fact (rescaling time and putting k, = 2 produces an agreement 
with our slightly different convention). The inductive definition of fk and the Sobolev 
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embedding theorem can now be used to show that, for £ > 2, 

(A.3) E2, = {-lY [ \d'hix)\^dx + Oiil + \\h\\^,^-if'^'), 

Jm. 

and hence for £ > 0, we have 

(A.4) \\d'hh^<c{i) 

where c{£) depends upon Sobolev norms of the initial data (f of at most order £. We now 
elaborate on how to obtain ()A.3|) . An inductive argument using ()A.2|) shows that for k >0, 
fk is of the form 

(A.5) h = hd'h + h J2 P(2j + l,A;-2j) 

j>l,2j<k 

where p{n, m) indicates a linear combination of terms of degree n and cumulative order m, 
or more precisely terms of the form 

(A.6) 9°i/i(9"2/i---9°"/i, ai + --- + a„ = m 

and h is either h or h. To prove ()A.3|) for £ > 2, one uses ()A.5|) for k = 2i and it only 
remains to verify that for any n > A and m < 2i — 2, 



J p{n,m)dx < \\h\\^t.i 



(A.7) 

We now show this. Note that in ()A.6|) . we may assume without loss of generality that 

Case 1. It follows that aj < £ — 2 for all j < n — 2 and a„_i < £ — 1. We 
estimate as: 



by Sobolev embedding. 

Case 2. a„ > £. In this case, we begin by integrating by parts to obtain 
(A.8) {-lY-^^+^ j d''"^^+\d'''h---d''-''h) d^-^hdx 

The Leibniz rule expansion is 

(A.9) =^c^9^i+"i/i---9^"-i+""-i/i 

where the sum is over (n — l)-tuples /i such that /ii + ■ ■ ■ + /in-i = — ^ + 1 and 
is some constant depending on /i. By adding the a and /i constraints, we obtain that 
(/Xi + q;i) + - ■ ■+ (/i„„i + a„_i) < i—1 and thus there is at most one index j^, (1 < < n—1) 
such that fij^ + aj^ = i — 1 and for all remaining j (1 < j < n — 1, j 7^ j*) we have 
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fij + Oij < i — 2. (If no such exists, take to be any fixed index 1 < < n — 1.) By 
substituting ()A.9|) into ()A.8|) . we estimate as 

f d'''hd'^^h---d''"hdx <( JJ ) < c||/i||^,_i 

again by Sobolev embedding. This concludes the proof of ()A.7|1 . thus ()A.3jl . and thus ()A.4jl . 

Using that the commutator [{x + itd),idt + ^d'^] = and some integration by parts 
manipulations, we have the pseudoconformal conservation law: 

f \{x + itd)h{x,t)\'^dx -t'^ j \h{x,t)\'^dx + j s j \h{x, s)\^dxds = j \x^{x)\^dx. 

J X J X Jo J X J X 

From this, (jA.4jl for £ = 0,1, and the Gagliardo-Nirenberg estimate ||/i||^4 < ||/i||^2 ||9/i||2,2, 
we have 

\\xh\\L2 < C{t) 

where c depends on ||xv9||^2, ||(/9||l2, and ||9(/?||j;,2. We want to show that more generally, for 
each G Z, /c > 0, we have 

(A.IO) \\x'^d'^h\\L2 <c{k){ty^^^ ioTa + (3 = k, a,/3>0, a,/3eZ 

Here c{k) is a constant depending on k and weighted Sobolev norms of the initial data ip 
(up to order 2k), and cr{k) is a positive exponent depending upon k. We are not concerned 
with obtaining the optimal value of (T{k); the mere fact that the bound in ()A.10|) is power- 
like in t, as opposed to exponential in t, suffices for our purposes. In our proof, both c{k) 
and cr{k) will be increasing with k, and will go to +oo a.s k ^ +oo. 

Let Aq = d and Ai = (x + itd). Note that both operators have the commutator property 



(A.ll) [A„izdt + ld^)]=0, j = 0,l. 

We first claim that for each A; > 0, there exists a constant c{k) > and an exponent 
a{k) > such that 

(A.12) ||A,,A,, ■ ■■A.MIl^ < c{k){ty^'^ for all j,, . . . j,. e {0, 1} . 

When we wish to consider a composition of the form Aj^Aj^ ■ ■ ■ Aj^, and do not care to 
report whether each operator in the composition is Aq or Ai, we will instead write the 
composition as A^. We prove ()A.12|) by induction on k . When k = 0, ()A.12|) is just the 
mass conservation law. Suppose that ()A.12|) holds for 0, . . . , /c — 1; we aim to prove it holds 
for k. The main ingredient (in addition to the inductive hypothesis) is ()A.4|) . Fix ji, . . .jk G 
{0, 1}, and apply the operator Aj^ ■ ■ ■ Aj^ to the equation, pair with —iAj^ ■ ■ ■ Aj^h, integrate 
in X, take twice the real part, and appeal to ()A.11|) to obtain 

(A.13) dt\\A,,---A,^h\\l2 = 2Rei I Aj,---A,Mh A,,---A,^hdx 
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Note that 

Ao F{h, h)h = dF{h, h)h + F{h, h) Aoh 

and 

Ai F{h, h)h = ttdF{h, h)h + F{h, h) Aih 
Both of these product rules take the form 

A \h\^h = g{t)dF{h, h)h + F(h, h) Ah 
where g{t) is a polynomial in t of degree < 1. Thus we see that 

A,,...A,Jh\'h=\h\%,---A,^h+ J2 9,it)d^'\h\' A^^h 

7i+72=fc 
72<fc-l 

where Q'yit) is a polynomial in t. Substituting into ()A.13|) . we obtain two terms: the first 
is zero since it is the real part of a purely imaginary number; the second is estimated by 
the Holder inequality to obtain: 

\dt\\Aj^---Aj^h\\l2\ <c{k){ty^^^ ( sup \\d^h\'^\\Lo.] ( sup ||A^'/i||i2 ) \\Aj^ ■ ■ ■ Aj^h\\L2 

\j<k-l J \j<k-l J 

By Sobolev embedding estimates, ()A.4|) . and the induction hypothesis, we have 
I 9,11 A,, ■ ■■A.Mll^ I < c(A;)(t)'^W||A,, ■ --A.MIl^ 

from which ()A.12|1 follows. 

Now to deduce (jA.lOjl from ()A.12|1 . we just note that since a; = Ai — it Aq, there are 
polynomials gj{t) such that the following relation holds: 



ie{o,i}°=+'3 



Let us now consider the application of (jA.lOj) to obtain in the Phase 3 analysis. 
We have xq + tv > v^^^ for t > ^2- If a; < 0, then 

Thus 

v''^'-'^\\utr{x,t)\\i^2^^ <\\{x-xo- t^;)'NLSo(t - t2)[t{v)sech]{x - xo - tv)\\Li 
= \\x''NLSo{t-h)[t{v)sech]{x)\\L2 
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by ()A.10|) . which gives the first estimate in ()3.9p . The second is obtained similarly. To 
obtain the third, we note that for t >t2, 

\utr{0,t)\^ = \NLSo{t - t2)[t{v)sech]{-xo - tv)\'^ 

dx\NLSo{t — t2)[t{v)sech.]{x — Xq — tv)\'^ dx 



and this can be estimated by 

||NLSo(t - t2Mv)sech]{x - Xq - tv)\\L2j\d,NLSo{t - t2)[t{v)sedi]{x - Xq - tv)\\Li^^ ■ 

Using ()A.10|) as before establishes 

^Mi-5)||a^NLSo(t - t2Mv)sech]{x - Xq - tv^^.^^ < c{k + j){t - ^2)^^^'+^'^ 

Replacing {c{k)c{k + l)y^'^ by c{k) and ^{^{k) + a{k + 1)) by a{k), we obtain the bound 
Finally, we note that the fourth bound in ()H.9j) is similar to the third. 

Appendix B. Free nonlinear evolution of asech 



This appendix is devoted to showing that for < a < 1 

NLSo((2a - l)sech((2a - !)•)) + OL^{t-^) 1/2 < a < 1 , 
CLoo(r^) 0<a<l/2. 



(B.l) NLSo(asech) 



A more precise understanding of error terms is possible thanks to advances in the study 
asymptotics for integrable nonlinear waves |3] , ^21 • Since we do not know an exact reference 
for ()B.1|1 we present a proof of this simpler asymptotic result. It is based on the now 
classical work on the inverse scattering method initiated for NLS by Zakharov-Shabat ^7j 
- see for discussion and references. For the reader's convenience, especially in view 

of different conventions used in different sources for our argument, we review all the needed 
aspects of the method. 

In the case of a = 1/2 we can use the result of ^21 to conclude that 
(B.2) NLSo(asech) = O {{log t/t)^) . 

A slightly inaccurate statement similar to (jB.ll) was given in ^Hj and the calculation of the 
scattering matrix in that paper was our starting point in obtaining (jB.ll) . 



B.l. Inverse scattering method. We present a quick review of this celebrated method. 
Thus, let us consider two operators acting on 5'(M; C^): 



L = -tJd, + tJQ, A = Jdl - ^ jg,. - JQd, - ^Q^J , 



where 



Q = Q{t,x) 



u{t,x) 



-u{t,x) 



u{t,»)ES(R), J 



-1 
1 
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which is checked by using 

JQJ-- 

It is now the case that 
(B.3) dtL = i[L,A] 



2 Qxx ~i~ 



■\u\ 




\u\ 



1 



idtu + -df.u+ \u\^u = 0, 



and, since we are solving NLS, we assume that these equivalent equations hold. 

We now consider scattering theory for the time dependent operator L. For that we 
introduce special solutions to Lip = \ip with prescribed asymptotic behaviour: 



(B.4) 



A) 



-ix\ 







-ix\ 







A) 

(^(x, A) 





^ix\ 



Ax\ 



X 



X 



+ 00 



-oo 



see instance [HI Sect. 1. 5]. Here for vector valued functions, ip =^ [</?2, —^it-, if = [v5i,v^2]*- 
Each pair of solutions forms a basis for the solution set and, for A G M, 

ip{x, A) = a{\)Tp{x, A) + h{X)ip{x, A) , 

(B.5) (^(x. A) = a{\)i){x, A) - h{\)i){x, A) , 

|a(A)r + |6(A)P = l. 

Another consequence comes from ()B.3|) . If L{t)il){t) = Xip{t) then we see that 

{L - X){tdtip - Aifj) =0 

and hence 

idtij{t) - Ai;{t) = ci{t)ij{t) + C2^{t) . 

Now we note that for u(t, •) G 5(]R), A ~ J9^, as |a;| — > oo, and the asymptotic behaviour 
()B.4|) gives ci(t) = A^, C2{t) = 0. More generally we conclude that 

idtip = iA + \'^)lP, idt^ = {A - X^)^ . 



(B.6) 



The solutions ip and ip have analytic extensions in A to the upper half plane and ip and <^ 
to the lower half plane. Same is true for a(A) and a(A) respectively. Except in very special 
cases (such as our potential a sechx) 6(A) does not have an analytic extension off the real 
axis. The reflection coefficient is defined as 

6(A) 



r(A) 



a(A)- 
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We assume that a(A) has at most one zero and that it can only he in Im A > 0, Re A = 0. 
That zero, Aq, corresponds to an LF' eigenfuction of L, and at A = Aq, the two solutions are 
proportional: 

(B.7) V5(x, Ao) = 7o^(a:, Ao) . 

The scattering data is given by the triple 

(B.8) u(t,a:) ^{r(A,t),Ao,7oW}. 



The evolution of the scattering data is easily obtained from (jR 



(B.9) 



r(A,t) 



r(A,0), 7o(t) 



7o(0) . 



In fact, we can use flB.5j) ()B.6|) to see 

{A — A^)v9 = idti^ = idtaip + aidtip + idtbijj + bidtijj 
= idtaifj + [idtb + 2\%)^l) + {A- \^)^ . 
Independence of ip and ip shows that, remarkably, 

9ta(A,t) = 0, dth{\t) = 2i\%\t) , 
which gives the first part of ()B.9|) . From ()B.7|) we see that. 



7o(A - A2)^(Ao) = idMXo) = idtio^iX,) + 7o(^ + A')V^(Ao) , 

so that, 

da^{t) = 2^270 (t) . 

That gives ()B.9j) . The justification of this formal calculation depends on u{t, 
and we refer to, for instance, Section 1.7] for a full proof. 



e S( 



B.2. The Riemann Hilbert problem. It is now universally acknowledged that the best 
way to obtain long time asymptotics for the inverse of ()B.8|) and ()B.9|) is by solving a 
Riemann-Hilbert problem ^ Chapter II]. To recall this method let us consider the 
following matrix valued function of A G C \ M, depending parametrically on x G M: 

[a(A)- V(a;, A)e'^^', ^{x, X)e-'^% Im A > 



(B.IO) 



dof 



[ip{x, A)e^^^, a(A)-V(a;, A)e-^^^], Im A < 
The properties of ip and (f (see for instance [3 Section 1.5]) imply that 



(B.ii; 



^(A) = / + 0(|Ar^), |A|^cx) 



where the decay rate may depend on x, uniformly in compact sets. From ()B.5|) we see that 

1 

I _f (_\)g-2i(Aa;+A2f) 
1 



^'(A + iO,x) 
^(A - iO,x) 



(x,A)e*^",^/^(x,A)e-'^"] 
(x,A)e^^^^/^(a;,A)e-^^^] 



A G 
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Hence, the boundary values of \I/(A) satisfy 

^(A + zO) = ^(A-zO)\4,t(A), AG 



(B.12) 



def 



x,t 



l + |r(A)|2 f(A)e-2^(^--+^'*) 
r(A)e 



2i{\x+\^t) 



If a(A) has no zeros in ImA > then the Riemann-Hilbert problem is to construct ^{X) 
satisfying (jB.llj) and ()B.12j) . Liouville's theorem readily shows that it is unique. If a(A) 
has a zero in ImA > 0, and in our presentation we allow at most one, Aq, we have to 
consider a Riemann-Hilbert problem in which \I/(A) is allowed to have singularities at Aq 
and Aq. The structure of that singularity can be seen in ()B.10|) : 



[B.13) 



ResA=Ao * = ^'(Ao 



ResA=Ao ^ = ^(Ao 





g2iAox^/^ 

e-2*^'o^7^ 



/ def 70 

7o 



a'(Ao) 



Since a(A) can be reconstructed from r(A) and Aq (see for instance Chapter I, (6.23)]; 
in our case it will be explicit) the Riemann-Hilbert problem in the case of one singularity 
is to find ^ which in addition to (jB.lljl and ()B.12jl satisfies ()B.18jl . 

A standard way to read off u{t, x) from \E'(A, x) follows from high frequency asymptotics 
of 'ip{x, A) (see for instance [T71 (18)]): 



(B.14) 

so that 
(B.15) 



il){x,\)e * 



1 





+ 



2i\ 



u(t, x) 
ir Ht,y)\^dy 



O 



u{x,t) = lim 2zA^i2(A,a;) 

A^oo 



We conclude this brief review by describing a reduction of the problem with prescribed 
singularities ()B.13|) to a problem with \l/ analytic in C \ M. To do that we follow |Sl Section 
II. 2] by considering a reformulation of the Riemann-Hilbert problem: 



G±(A)^i = *(A) 



1 

a±(A)^i 



a+(A) =^ a(A) , a_(A) =^a(A), ±ImA>0 



G+(A + iO)G_(A-20) =G(A), G(A) 



def 



1 -6(A) 
-6(A) 1 



The operators G±{\) are now analytic in ImA± > (in fact, ^^^(A) = G^{X)) but their 
ranks drop precisely at A = Aq and Aq respectively. The condition ()B.13|) becomes 



(B.16) 



ImG_i_(Ao) = span^ 



7oe 



1 

2iXox 



kerG'_(Ao) = span^ 



-7oe" 
1 



-2i\QX 
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We now look for B{X), analytic in IniA > 0, with B{X) ^ analytic in ImA < 0, -B(A) = 
J + C(1/|A|), and such that 

G+{X) = G+{X)B{Xy\ G_(A) = 5(A)G+(A), ±ImA>0, 

are nonsingular matrices. We note this requires -B(A)~^ to have a pole at Aq and B{X), at 
Aq. That is natural since we are adding to the ranks of G±{X). The condition ()B.16jl mean 
that 
(B.17) 

1 



Imi?(Ao) = spanj;.^^ (Aq 



keri?(Ao) = spancG'_(Ao) 



-7oe 
1 



-2iXox 



This determines -B(A) uniquely as a Blaschke-Potapov factor: 

, Ao — Ao 



BiX) 



(B.18) 



P 



1 



W 

p 1 



A-Ao 

(3{x) 



-P, P* 



P. 



GL"^(Ao,x)7oe-2^^o- + G!_''^(Ao,x) 



(12), 



GL'')(Ao,x)7oe-2^^o- + GL^^^(Ao,x)' 
see Chapter II, (2.17)-(2.27)]. 

Hence to solve the Riemann-Hilbert problem (jB.12jl . (jB.13|) we first solve the problem 

(B.19) $(A + zO) = $(A-zO)\4,i(A), ^ = / + C(1/|A|) , $ analytic in C \ M . 
We then define 

" 1 



^(22), 



G_(A) = ^(A) 



a(A) 



G+(A) =G_(A)*, ImA<0^ 



from which, using ()B.18|1 . we construct B{X). Then 
(B.20) *(A) = fi(A)-^^(A), 

and we can finally use ()B.15|) to obtain u{t,x). In particular the long time behaviour of 
u{t,x) is determined by the longtime behaviour of ^(A). 

B.3. Manakov ansatz. The basic structure of the long time behaviour of \E' solving (jB.19j) 
can be obtained from the Manakov ansatz for the solution of ()B.19jl - see [H Section 2] and 
references given there. To describe it we define 



m±(A,x,t)^=^^ 



±1 



where r+(C) = r(C), r_(C) = r(C), and 



6{X,x,t) = exp 



log(l + |r(C)p- 
C-A 



dC. 



dC 
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solves a scalar Riemann-Hilbert problem, 

' \ 5(C-«0), C>-x/2t, 

see jSl Proposition 2.12] for a detailed list of properties of 5{z) (stated in the defocusing 
case log(l — l^'(C)P))- The Manakov ansatz is then given by 



(B.21) 



*(A,x,t) 



def 



1 m+(A,x, t) 
m„(A,x, t) 1 



(5(A,x,t) 1 
1 5(A,x,t)-i 



To see the properties of ^ we use the following elementary lemma: 

Lemma B.l. Suppose that f e C°°((0,M)), /' G L\[0,oo), x^/^'^ G L°°([l,oo)), for all k 
and I. Then, as X ±oo, 



fiy)e 



27ri J Q y — X — iO 
Using this lemma one checks that 



f(^rj.y\x^ + C(|A|-5|x|-i) x > arw^ A > 
(9(|A|~5 otherwise. 



from which it follows as in jH Section 2] that 

$(A) = (J + 0(lM))$(A). 

If we defined ^(tjx) by putting \& in ()B.15jl this shows that u{t,x) = 0{l/\/i), and in fact 
a more precise statement can be obtained by using the second component in ()B.14|) . 

We can now use ()B.18|) to obtain an approximation to -B(A) as t ^ oo: 
B{X) = 5(A) + 0{t~'/^\X\^^) , B{X) = 1- ^!LZ^p , 



A-Ar 



where P is as in ()B.18|) with f3 replaced by 
/3 = 7o exp(— 2Aox) exp 



1 

iri 



° log(l + |r(C)p 

oo C ~ -^0 



-dC> . 



Hence to obtain a long time approximation for a solution we apply the procedure of jSl 
Section II. 5] to Bo[X) since that corresponds to using (IrTbI) with 

^(A) = B{xy^^{x) = B{x)-\i + o(ri/'))$(A) = Bo{xy\i + c(r^/'|Ar^)) . 

This gives 
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Lemma B.2. Suppose that u{0, •) G S{M.) and that the scattering data ()B.8|) for u{0, •) is 
given by r{X), 70 and Aq, ImA > 0. Then 



(B.22) u{x,t) = e^-^oNLSo (e2^^"^°'2ImAosech(2ImAo(« -xo))) +CL-(r 



l/2^ 



where 



= 21^ 1^°" - I" - log(2ImAo) - j ^ (C-ReAoP + Ag ^^ 

and 

1 /"^ log(l + |r(C)P) 
ifQ = arg7o - arga'(Ao) + - / 77 — p x no , \2 (C - ReAo)rfC- 

i_oo (C - ReAo)^ + A^ 

We state this important result as a lemma to stress the fact that a better error estimates 
seem available if more advanced methods are used. 



B.4. Scattering matrix. We now apply Lemma IB .21 to obtain ()B.1|) . For that we need 
to find the scattering data ()B.8|) for u{0, x) = asechx. That is done by a well known 
computation ^HIj^I Sect. 3. 4] which reappears in many scattering theories, from the free 
S'-matrix in automorphic scattering, to Eckhardt barriers in quantum chemistry. We quote 
the results: 

r(i-zA) . sinvra f _ K>^) 

"^^^-r(| + «-zA)r(i-«-zA)' ^^^'-'Z^sh^' 

We note that in this special case b and r are meromorphic in C (with infinitely many 
"nonphysical" poles). Also, 

Ao = 2a — 1 if 1/2 < a < 1 , and 70 = &(Ao) = i ■ 

We need to compute xq and ipo. In general when u{0,x) is even then xq = by symmetry 
considerations. Here we see it by using Chapter II, (2.6)] which shows that 

ImAo Z"" log(l + |r(C)r 



TT C^ + ImA; 



2 



-dC. 

Thus the formula in Lemma IB. 21 results in xq = 0. To compute (fo we need to find the 
following integral 

/ sin^ na \ C , ^ , 
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